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The twisted reduced model of large TV QCD with two adjoint Wilson fermions is studied 
numerically using the Hybrid Monte Carlo method. This is the one-site model, whose 
large N limit (large volume limit) is expected to be conformal or nearly conformal. The 
string tension calculated at A'^=289 approaches zero as we decrease quark mass and the 
preliminary value of the mass anomalous dimension 7* is close to one if we assume that 
the theory is governed by an infrared fixed point. We also discuss the twisted reduced 
model with single adjoint Wilson fermion. The string tension remains finite as the quark 
mass decreases to zero, supporting that this is the confining theory. 
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1. Introduction 

The standard model of the elementary particle has set the foundation to SU(N) 
gauge theories. Although SU(N) gauge theories generally have very complicated 
C^ ■ structures, significant simplifications could occur in the large N limit. In fact, Eguchi 

and Kawai considered the space-time reduced model of the lattice gauge theory hav- 
ing only one space-time point. ^ This model is now called the Eguchi-Kawai model 
(EK- model). The EK model possesses the Z(N) symmetry. If this symmetry is not 
spontaneously broken, the lattice gauge theory and the EK model are equivalent 
in the large N limit. This means that space-time degrees of freedom of the lattice 
gauge theory could be encoded into those of internal SU(N) group. Soon after the 
proposal of Eguchi and Kawai, however, the Z(N) symmetry was shown to be spon- 
taneously broken in the weak coupling limit, thus invalidating the EK model. ^ To 
circumvent this difficulty, the present authors have proposed to consider the one-site 
model with twisted boundary conditions (the twisted Eguch-Kawai (TEK) model), 
thus avoiding the Z(N) symmetry being spontaneously broken.^' Recently, gauge 
theories coupled with fermions in the adjoint representation received much attention 
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in the large N limit in connection, for example, to the AdS/CFT correspondence. 
The idea of the twisted space-time reduction can also be applied to these models. 
The purpose of this talk is to present our recent progress in this direction. In the 
next section, we calculate the continuum string tension of the SU(N) pure gauge 
theory in the large N limit using the TEK model. In sect. 3, large N QCD with two 
adjoint Wilson fermions is studied with the twisted space-time reduction, and sect. 
4 is devoted to the large N QCD with single adjoint Wilson fermion. Conclusions 
are given in sect. 5. 

2. The TEK model 

We consider the SU(N) group with N ^ L"^, L being some positive integer. Then 
the action of the TEK model is obtained from the action of the SU(N) lattice gauge 
theory 

4 

S^-bNY, E Tr [t/„,^t/„+^,,C/l+,^^C/^,,] (1) 

n ^5^jy=l 

by neglecting the space-time dependence of the link variable Un,fi ^- Uf^ and multi- 
plying the twist tensor Z^^ in front of the plaquette action as"^ 

4 

STEK = -bN J2 Tr [Z^.C/^t/.t/^C/t] (2) 

As a result, there remain only four SU(N) matrices U^. The inverse 't Hooft coupling 
is labelled b = l/g^N . The elements of the twist tensor belong to Z(L) and the 
explicit form is given by 

Z^^ = exp f k— j , Z^^ = Z*^, 11 > V (3) 

The action Stek is invariant under the Z(L) transformation [/^ — > zU^, z G Z(L). If 
this symmetry is not spontaneously broken, the TEK model and the corresponding 
lattice gauge theory are equivalent in the large A'' limit, fc and L should be coprime, 
and a general prescription for choosing fc and L to preserve the Z(L) symmetry has 
been given in ref. [4]. In this talk, we choose L=29 and fc=ll for the TEK model. ^'^ 
For finite L, the TEK model is closely related to the lattice gauge theory on a finite 
space-time volume L**, and we are able to calculate Wilson loops up to size 14 x 14 
for L=29. In the TEK model, the Wilson loop W{R, T) with size RxT is obtained 
from four link variables C/^ as 

W{R, T) = Z^^ < U^U^U^^Ur > ■ (4) 

Then the string tension a is extracted from the large distance behavior of the Creutz 
ratio 
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Fig. 1. A'^ dependence of the continuum string tension. 
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as 



x(i?',i?')^. + |j + ^ 



(5) 



(6) 



with half-integer R' and T' . We calculate the continuum string tension by extrap- 
olating the TEK data at six values of b. For comparison, we also calculate the con- 
tinuum string tension of the usual SU(N) lattice gauge theory with iV = 3, 4, 5, 6, 8 
on a y = 32'' lattice. Results are summarized in Fig. 1. The x point is Ajj-g/^/a 
for the TEK model with 7V=841. The -I- symbols are results for the lattice gauge 
theory, and dotted line is a linear fit of these data as a function of l/N"^ . The large 
N extrapolated value agrees remarkably well with the result of the TEK model, 
demonstrating the correctness of the twisted space-time reduction idea. 

3. Large N QCD with two adjoint fermions 

SU(N) gauge theories with two adjoint fermions are thought to be conformal or 
nearly conformal for any value of iV, since the first two coefficients of the beta 
function expressed in term of 't Hooft coupling are independent of N , namely, 
6o = (4iV/ - ll)/247r2 and fei = (167V/ - 17)/1927r4. For iV/ = 2, &o < and 
bi > 0, then we naturally expect that there is a infrared fixed point for finite value 
of 't Hooft coupling. In fact, for N=2 (minimal walking technicolor), there are 
now many lattice simulations indicating that the theory is conformal at vanishing 
fermion mass.^ The mass anomalous dimension, however, is shown to be rather 
small 7* ^ 0.3, to explain the large values of the observed quark masses, within the 
walking technicolor scenario. 
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Siring tension and lowest eigenvalue of Q, N,=2, SU(289), k=5, b=0.35 
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Fig. 2. String tension and lowest eigenvalue of Q^ at 6=0.35 with A'^ = 2. 

We can also consider the twisted space-time reduced model of large TV QCD 
with two adjoint fermions. |fj The action reads 



4 Nj ^ 

S = -bN ^ Tt[Z^,U^U,uIuI] +Y,^jDw^j 
with Nf=2 and the Wilson-Dirac operator Dw given by 



(7) 



z?w = 1 - /. ^ [(1 ~ j,)uf^ + (1 + i.Wl^'q 



M=l 



(8) 



"^j is the fermion matrix in the color (N, N) representation. Thus the link variable 
in the adjoint representation C/°* actually acts on ^j as C/°*^j ~ U^'^jUJ^. The 
hopping parameter of the fermion field n is related to the bare fermion mass as 
TO/ = (1/2)(1/k— 1/kc) with Kc the critical hopping parameter. 

We calculate the string tension at two values of 5, 0.35 and 0.36, and for various 
values of k. We are using N — 289 — YI^ and fc = 5, which should compare with 
ordinary lattice theory with V — 17^}^ For Nf = 2, we can use the standard hybrid 
Monte Carlo method. Simulations have been done on Hitachi SR16000 supercom- 
puter at KEK having peek speed of 980 GFlops per node. Thanks to the Hitachi 
system engineers, our code is highly optimized for the SR16000. The sustained speed 
of our code being 600 GFlops. 



''For previous works of the space-time reduced model of large A^ QCD with adjoint fermions, see 
ref. [8-10]. 
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String tension and lowest eigenvalue of Q^, N,=2, SU(289), k=5. b=0.36 
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Fig. 3. String tension and lowest eigenvalue of Q^ at 6=0.36 with Nf = 2. 

In Fig. 2, we show the string tension cr as a function of k at 6=0.35 with + 
symbols. The value of a at k=0 is obtained from the TEK model without fermions. 
We clearly see that as k increases the string tension rapidly decreases and seems 
to vanish at k '-^ 0.17. So far, we have not calculated any hadronic spectrum. It 
is quite straightforward, however, to calculate the lowest eigenvalue of the positive 
hermitian Wilson-Dirac operator Q^ = {Dwj5)'^, which should be related to the 
physical fermion mass square. In Fig. 2, we also show the lowest eigenvalue of Q^ as 
X symbols. The string tension and the lowest eigenvalue of Q^ vanish simultaneously, 
which strongly supports that the string tension is zero at critical Kc where the 
fermion mass vanishes. 

We can fit both the string tension a and the lowest eigenvalue of Q^ with the 
same fitting form a(l/K— 1/ Kc)^. Requiring that both quantities vanish at the same 
critical Kc, we have Kc — 0.1694(7) and, for the string tension, b = 1.10(4). If the 
massless theory is governed by an infrared fixed point, b is related to the mass 
anomalous dimension 7* as 6 = 2/(1 + 7*), then we have 7, = 0.81(8). 

We have repeated the analysis at &=0.36, a little bit closer to the continuum 
limit. Results are shown in Fig. 3. It seems that the data at k = 0.16 are largely 
affected by the finite -size effects. Excluding the data at k = 0.16 from our analysis, 
we have Kc = 0.168(1), b — 1.33(4), and thus 7* — 1.17(21). Although our data are 
still preliminary, large N QCD with two adjoint fermions seems to be a conformal 
field theory with large mass anomalous dimension 7* ^ 1. 



4. Large N QCD with single adjoint fermion 

In the large N limit, Nf — 1 adjoint fermion is equivalent to the Nf — 2 fun- 
damental fermions in the rank two anti-symmetric representation (the orientifold 
equivalence*). For A^ = 3, the latter theory is just the two flavor QCD. This means 
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String tension and lowest eigenvalue of Q , Npl , SU(2S9), k=5, b=0.35 
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Fig. 4. String tension and lowest eigenvalue of Q^ at 6=0.35 with Nf = 1. 



String tension and lowest eigenvalue of Q, N,=1, SLJ(289), k=5, b=0.36 
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Fig. 5. String tension and lowest eigenvalue of Q^ at 6=0.36 with Nf = 1. 



that the large N QCD with single adjoint fermion corresponds to the Corrigan- 
Ramond large N limit^^ of realistic QCD, and thus it is quite interesting to analyze 
this theory from the phenomenological point of view. 

We use the rational hybrid Monte Carlo method to simulate the Nf = 1 theory. 
In Fig. 4, we plot the string tension cr as a function of k at 6=0.35 with + symbols. 
Also shown is the lowest eigenvalue of Q^ with x symbols. The lowest eigenvalue 
of Q^ seems to vanish around k ~ 0.165, while the string tension seems to remain 
finite at this point, indicating that the large N QCD with single adjoint fermion is 
a confining theory. This is quite natural since our theory should be related to the 
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two flavor A^ = 3 QCD. We found ttiat, for k > 0.16, we cannot make simulations 
since the CG iteration does not converge, which is also the common phenomena 
encountered in the simulation of QCD. 

In Fig. 5, we show the results at 6=0.36. The absolute values of the string tension 
are smaller than those of 6 = 0.35 since we are approaching to the continuum limit. 
The string tension, however, remains finite around k ^ 0.165 as in Fig. 4. We also 
confirm that the CG iteration docs not converge for k > 0.16. Hence, both Figs. 4 
and 5 indicate that the continuum theory of the large N QCD with single adjoint 
fermion is a confining theory. 

5. Conclusion 

We have demonstrated that the twisted reduction works quite well for the descrip- 
tion of large N QCD. For pure gauge theory, we have succeeded in calculating the 
continuum string tension using TEK model. Although our calculations with dy- 
namical fermions are still preliminary, we have also shown that large N QCD with 
two adjoint fermions is a conformal theory with large mass anomalous dimension 
7, ^ 1, and that the theory with single adjoint fermion is a confining theory. Since 
the rational hybrid Monte Carlo method we used for the simulation with one fla- 
vor can directly applicable to any number of flavor, we are planning to extend our 
analysis to Nf = 1/2 super-symmetric theory and also Nf = 3/2 and 5/2 theories 
if they are conformal or confining theories. 
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